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Abstract 

In this work we investigate the issue of gravity and fermion localization and resonances in (4, 1)- 
branes constructed with one scalar field coupled with gravity in deformed models. Such models 
give solutions for the scalar field that is the usual kink solution in the extra dimension for a 
parameter p = 1 and deformations with a two- kink profile for odd p > 1. Gravity is localized 
and resonant modes are found for small values of p. The coupling between the scalar field and 
spinors is a necessary condition for fermions to be localized on such branes. After performing a 
chiral decomposition of the five-dimensional spinor we found resonances with both chiralities for 
all odd p's. The correspondence between the spectra for left and right chirality is guaranteed and 
Dirac fermions are realized on the brane. The increasing of p characterizes the formation of branes 
with internal structure that turns the gravitational interaction more effective for fermions aside 
the brane, increasing their lifetime. The influence of the internal structure of the branes and the 
presence of resonances for gravity and fermionic modes is addressed. 

PACS numbers: ll.10.Kk, 03.50.-z, 04.50.-h, 11.27.+d 
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I. INTRODUCTION 



Brane structures where initially introduced as domain walls embedded in extra dimensions 



naa 



In a sense thick branes 



are a natural construction since the solutions 
can be dynamically found. More recently, interesting models have been prop osed for such 
dynamical branes with rich structures, constructed with one [uj] or more ll[ 12, 13] scalar 
fields (see also the review given by Ref. [lij]). The issue of localization of several fields and 
resonances in such branes is an interesting subject, as their investigation can guide us to 
which kind of brane structure is more acceptable phenomenologically. In this way, models 



which extend the Randall-Sundrum type II scenario 
under the aspect of gravity localization jl^l . 



17 



18 



19 



5| were constructed and considered 



20]. 



Reviewing works about fermions in extra dimension models, we observed that in type I 
Randall-Sundrum (RS) scenario 21(, spin 1/2 and spin 3/2 fermions, are located only in 
the negative tension brane. However, they do not satisfy the localization conditions in the 
type II RS model. In these two cited models, the localization of zero mode chiral fermions is 
only obtained due to the introduction of a generalized Yukawa coupling. This is analogous 
to what happens to domain walls in the absence of gravity 22j. In the brane case the 
coupling is introduced by an interaction of fermions with a bulk scalar field. The scalar field 
solution can be kink like, for instance. Thus, the results with the RS model indicate that the 
thick brane scenarios as the most suitable for locating fermions. Thick brane models can be 
obtained from scalar fields with kink-like solutions. Therefore, the scalar field that couples 
with fermions in the Yukawa coupling will be the scalar field which the brane is made of. 



Quite recently, Liu et al. 



231 ] have analyzed the issue of fermion localization on a pure 



geometrical thick brane. Since then, the same group of authors have published a series of 
papers on fermion localization on several types of thick branes 



24 



25 



26 



271 ] . Also, some 



of the authors of the present work were the first to consider fermionic resonances in branes 
with internal structure 

As is well known, the kind of structure of the considered brane is very important and 
will produce implications concer ning the methods of field localization. In the seminal works 
of Bazeia and collaborators [29j, |30J a class of topological defect solutions was constructed 
starting from a specific deformation of the <f) A potential. These new solutions may be used 
to mimic new brane- worlds containing internal structures 3j|. Such internal structures 
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have implications in the density of matter-energy along the extra dimensions 11[ and this 
produces a space-time background whose curvature has a splitting, as we will show in this 
work, if compared to the usual models. Some characteristics of such model were considered 

n 

in phase transitions in warped geometries [8J. 

Considering a brane world scenario like a 4D domain wall immersed in a 5D space-time, 
we were able to find a localized left chiral zero mode. Therefore, if we consider a deformed 
brane, we obtain new results related to the localization of spin 1/2 fields. Also, such class of 
models where already studied with respect to gravity localization [10|, where it was found 
that zero-modes for KK gravitons exist, the solutions are stable and tachyons are forbidden. 

The main goal here is to study the behavior of fermions and gravity in a membrane with 
internal structure generated by a deformation procedure. We modify the so-called two-kink 
solutions that can be obtained after a deformation procedure of a potential from a scalar 
field y. We analyze fermion localization on such branes, considering such deformations 
that suggest the existence of an internal structure. As we will see, the deformations will be 
very important for localization and normalization of fermionic fields. Using a well known 



resonance detecting method 26|, |28|, |3j|, we analyze the massive modes arising from the 



dimensional reduction functions. Modes with large amplitudes near the brane can reveal us 
details about the coupling of modes with the brane. 

This paper is divided as follows. In the next section we review the deformed brane model 
with one scalar field, where we analyze the bulk scalar curvature, relating the results to 
the the brane internal structure. Sec. Ill is devoted to the analysis of gravity localization 
and resonances with KK gravitons. Each particular brane solution is considered as a fixed 
background that are not perturbed significantly by the presence of fermions. With this 
approximation, in Sec. IV we perform a chiral decomposition and study the presence of zero 
modes with left chirality. Massive modes and their normalization are subject of investigation 
in Sec. V. Fermionic resonances and the realization of Dirac fermions on the brane are 
considered in Sec. VI. The main conclusions are presented in Sec. VII. 
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II. BRANE SETUP 



We start with the action describing one scalar field minimally coupled with gravity in 
five dimensions 

S = J d b x^G[2M 3 R - ^(<90) 2 - V (</>)], (1) 

where M is the Planck constant in D = 5 dimensions and R is the scalar curvature. For some 
classes the potential V((f>), it is possible to obtain kink solutions for the field depending 
only on the extra dimension. As an Ansatz for the metric we consider an extension for 
the Randall-Sundrum metric, where the bulk spacetime is asymptotically AdS$, with a 
Minkowski brane, 

ds 2 = e^Mri^dafda? + dy 2 (2) 

The warp factor depends on the metric function A(y), where y is the extra dimension. The 
tensor r]^ is the Minkowski metric and the indices fi and v vary from to 3. For this 
background we find the following equations of motion: 

#" + 4A V =g, (3) 



1 

2« 



^ 2 - V{4>) = 24M 3 (A') 2 , (4) 



and 

1 



(0') 2 + V{<j>) = -12M 3 A" - 2AM 3 (A') 2 . (5) 



Here prime means derivative with respect to the extra dimension. 
In the presence of gravity, defining the potential as 

1 fdW\ 2 8M 3 



W) = i(^r) ~^w 2 , (6) 



it is possible to find first-order equations 

dW 



(7) 



W = -3A'(y), (8) 

whose solutions are also solutions from the equations of motion. Here W(<j>) is the superpo- 
tential. This formalism was initially introduced in the study of non-supersymmetric domain 



walls in various dimensions 
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For bounce-like solutions, the field tends to different values when y — > ±00. Such 



so 

331 



utions can be attained by a double-well potential. In this way, guided by refs. 



10 



291 



341 ]. we chosen the superpotential, 



p 



2p-l 
) P 



P 



2p+l 
) P 



(9) 



2p - 1 T 2p + 1 

where the parameter p is an odd integer. The chosen form for W p was constructed after 
deforming the A0 4 model. This choice allows us to obtain well-defined models when p = 
1, 3, 5, where for p — 1 we get the standard 4 potential. For p — 3, 5, 7, the potential 
V p presents a minimum at = and two more minima at ±1. Eq. j7j) can be easily solved 
giving the so called two-kink solutions 



,(y) = tanh p 



(10) 



From Eq. (jSj), one find explicitly the solution for A p (y) as [lOj 



= -i^tanh 2 * - ? f-^- 

pVy; 3 2j9 + i V^y 3 V 2 p — 1 



p- 



2p 



(11) 



111 



ooahf^l -V-tanh 2n 




FIG. 1: Plots of the solution of the curvature invariant R{y) for p = 1 on the left. On the right for 
p = 3 (doted line), p = 5 (dashed line) and p = 7 (solid line). 



Note that the exponential warp factor constructed with this function is localized around 
the membrane and for large y it approximates the Randall-Sundrum solution 151 ] . The 
spacetime now has no singularities as we get a smooth warp factor (because of this, the 
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model is more realistic) 35[]. This can be seen by calculating the curvature invariants for 



this geometry. For instance, the Ricci scalar is 

R = -[8A; + 20(A> p )% (12) 

Figs. [T^, and[Tb shows the Ricci scalar for p = 1,3, 5, 7. Note that the Ricci scalar is finite at 
all points in the bulk, and singularities are absent. Far from the brane R tends to a negative 
constant, characterizing the AdS§ limit from the bulk. Note that the higher is the parameter 
p, the lower is the constant, interpreted as the inverse of the AdS scale. Fig. [Th is for p — 1 
and shows that the scalar curvature for a non-deformed model presents a maximum at the 
brane center y — 0, whereas for higher values of p this maximum is splitted into two smaller 
peaks that decrease with the increasing of p. The presence of regions with positive Ricci 
scalar can in principle be connected to the capability to trap massive states near to the 
brane, as we will investigate in the following sections. Also note that for even larger values 
of p we see that the tendency is the Ricci scalar to approach to zero near to the brane center. 
In this way we expect the presence of gravity resonances to be more pronounced for lower 
values of p. 

III. GRAVITY LOCALIZATION AND RESONANCES 

The issue of gravity localization in this class of models was considered by one of us in 



Ref. [lOj. As already noted in [7], plane wave solutions of Schro dinger- like equations in the 
transverse-traceless sector of metric perturbations can present solutions as resonant modes. 
Such structures were obtained in {(j, \3(^ when studying gravity localization. Here we review 
the treatment of metric fluctuations and investigate numerically the presence of resonances 
with the more refined method. The stability analysis is performed after perturbing the 
metric as follows 

ds 2 = e 2Afo) (r^ + h ilv )dx ll dx v - dy 2 . (13) 

Here = hfj, v (x,y) are small perturbations. In the transverse traceless gauge the pertur- 
bations turn to h^ u , and the metric and scalar field fluctuations decouple, resulting in the 
equation 

h»+AA'h'=e- 2A nh lu ,. (14) 
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FIG. 2: (a) (upper left) Plots of the potentials U\(z), U%(z), an&U^z). Normalized \H^ U (0)\ 2 , as 
a function of m, for m > 0, showing the resonance peaks for (b) p = 1 (upper right), p = 3 (lower 
left) and p = 5 (lower right). 



Here □ stands for the 4- dimensional D'Alembertian. The extra dimension y is turned into 
a new coordinate z, defined by 

(15) 



dy 



e~ Ap , 



which makes the metric conformally flat. Also, with 



(16) 



the equation (j!4p for the metric fluctuations assume the form of a Schrodinger-like equation 

d 2 H 

a i t . . . . , .-, . . (17) 



where the potential is given by 



U p {z) = \AZ{z) + \j%{z). 



(18) 



Is was already shown [lCj that the Hamiltonian is positive definite and tachyonic modes 
are absent and that it is possible to attain an explicit expression for the non-normalized zero- 
modes, responsible for gravity localization. Fig. [2^ shows that the Schrodinger potentials 



for gravity fluctuations have the form of volcano potentials. This inspired us to investigate 



the possibility of resonances with the Numerov method [18] . identifying resonances as a 
peak in the normalized squared wavef unctions \H^ V {Q)\ 2 at the brane center. We considered 
normalization in a box with ends at ±z max , far enough for the inverse square law U p (z) ~ 
a p (a p + l)/z 2 to be achieved. As gravity localization can be determined by the far region 
of the potential Q, Q, the plot of z 2 U p (z) for -200 < z < 200 gives e*i = 1.483. This 
characterizes gravity localization where U grav (r), the gravitational potential between two 
unit masses distant r one from the other, reproduces the Newtonian limit for large distances. 
However, for short distances there is a l/r 2ap correction due to the small massive modes. For 
p = 1 this gives l/r 2 966 , close to the Randall- Sundrum 1/r 3 correction for the Newtonian 
potential. For larger values of p, one needs larger values of z max to achieve the 1/z 2 region 
for U p . 

For studying resonances, we do not need to pursue to large values of z max , since we are 
interested in the effect where the wavefunction behavior changes abruptly for a particular 
mass. For our purposes we considered sufficient to consider a normalization procedure with 
Zmax = 100. We found a clear peak for p = 1, as showed in Fig. [2b. We noted that with the 
increasing of p (Fig. [2b and[2bl), the resonance peaks become broader, showing that branes 
with smaller values of p are more effective in trapping graviton KK modes. This can also be 
noted with zero modes, where the modes decay as 1/ yfz for z > z p and z p grow with p (see 
Ref. jl0|). For p = 5 the peak thickness Am is too large for characterizing a resonance. This 
was expected since the maxima of the Schrodinger potential, firstly pronounced for p = 1 is 
reduced considerably for larger values of p (see Fig. EK) • 



IV. FERMIONIC ZERO-MODE 



Now we consider the action for a fermion coupled with gravity, in the background given 
by the brane solution 4> p (y) from the previous section. The action is 

S = J dx 5 ^r M D M V - (19) 

where / is the 5-dimensional Yukawa coupling. 

Following the literature, firstly we change the variable y to z according to Eq. (|T5|) . In 
the new set of variables the metric is conformally plane, and the gamma matrices can be 
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rewritten as T M = e a i>j^ , r 5 = e a p^ 5 . The covariant derivatives are 

fl M = 3m + n^Vr 6 , ^5 = 5 5 . (20) 
The former expressions allow us to write the equation of motion as 

[7^ + l\d z + 2d z A p ) + f<P p e A ^(x, z) = 0. (21) 
Now we perform a chiral decomposition of the 5-dimensional spinor \1/ as 

2;) = ^[fe(i)a Ln + VijnO^aRnW]- (22) 

n 

The massive modes from the spinor \I/ living on the brane must connect both chiralities, 
satisfying the equations 

Yd^Lnix) = mip Rn {x) , Yd^Rn(x) = mip Ln (x). (23) 

From the relations "f 5 ipLn( x ) = —ipLni.%) and ^ipRnix) = tpRn(x), we find two coupled 
equations for ai n (z) and a Rn (z): 

[d z + 2d z A p + f<j) p e Ap }a Ln (z) = ma Rn (z) } (24) 

[d z + 2d z A p - f<p p e A v]a Rn (z) = -ma Ln (z). (25) 

Now we investigate the possibility of localized fermionic zero massive modes. For m = 0, 
Eq. (1241) for a>L n (z), reduces to 

2A' p a Ln (z) + a' Ln (z) + fe A ^ p a Ln (z) = 0. (26) 

If we turn back to the variable y we can write 

2A' p a Ln (y) + <x r Ln(y) + f^Ln(y) = 0. (27) 

with solution 

a Ln (y) = e~^ dyV<t>p+2AM . (28) 

In this solution we clearly note the contribution of the internal structure from the membrane. 
Remembering that 4> p and A p depend on odd integer numbers, we will see that the value of 
p will be determinant in order to obtain a finite solution. With the explicit expressions for 



A(y) and 4> p {y) in Eq. (|28|) . and following Ref. [35j, we where able to find the relation 
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between the coupling constant / and the parameter p for aL n (y) to be finite. 

As an example, we consider the particular case where p = 3 in order to determinate the 
effects of the deformations on the issue of the localization of the solution. In this case, Eq. 
can be written as 



24 



aLn ^ = e -y S ech(^+^tanhqr[-6-3tanh^)Hmanh(^] msh 35 3 / _ ^ 

The exponential part of the solution above tends to a constant value in regions far from 
y = for any value of /. Therefore, the exponent |g — 3/ determines the form of the 
solution when y — > ±00. In order to obtain a finite solution a following condition must be 
obeyed 

24 

35 " 3/ < °- (31) 

Analyzing the solution for A p (y) in ( fill , we note that the factor 24/35 in the solution for 
p = 3 results from the term, 

U p 2 p 2 \ (32) 



3 \2p- 1 2p+ 1 # 

Simplifying the term above, we obtain the general relationship from Eq. fl29l) between / e p 
that must be satisfied for all solutions for ai n {y) to be finite. It is important to point out 
that the relation from Eq. (T291) between / e p results from the coupling of the fermion \1/ with 
the particular kind of membrane solution, introduced in the action as f^f(f) p ^f. Fig. (jHfe.) 
shows the solution amiy) for p = 1,3,5 for y > 0. Note from the plots that the solutions 
asymptote to zero for large y, as required for localized solutions. Also note that small values 
of p shows higher peaks, characterizing better localization. 

For a Rn (y), Eq. ff25l) for m — 0, after turning to the variable y, reduces to 

2A' p a Rn (y) + a' Rn (y) - f(j) p a Rn (y) = 0, (33) 

with solution 

aBn (y)=er«V*>W-™,W)) m (34) 

From the solutions A(y), 4> p (y) we conclude that there is no localized right chiral zero-mode 
for / satisfying Eq. f[2T?j) . This can be seen explicitly in Fig. (j3b) where we plot solutions 
from Eq. ( |34"1) for p = 1,3,5. Note from the figure that the p — 1 case diverges for lower 
values of y in comparison with larger values of p. At least for zero modes we conclude that 
lower values of p favor localization of left chiral modes and disfavor localization of right 
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FIG. 3: (a) Plots of zero-mode aL n (left) and a Rn (right) with / = 1 for p = 1 (dashed line), p = 3 
(doted line) and p = 5 (solid line). 

chiral modes. We expect this to be also valid for massive modes. Further we will confront 
our findings for the zero-modes with the general calculations for the massive modes. 

Now to complete our analysis for chiral zero-modes we must verify the normalizability of 
our solutions. Decomposing the y-dependent part from the action for free fermions given by 
the first term from Eq. (|T9l) we get 



' free 



/r+oo 
d 4 x / dyy/g*(x,y)r A D A *(x,y) 
J — oo 

dye 3Ap{y) \a(y)\ 2 / d A xi>(x)YdM x )- 



(35) 



This expression shows that only left zero-mode chiral solutions are normalizable, as lead to 
a finite integral in y 351 ] . 



V. FERMIONIC MASSIVE MODES 



Now to complete our investigation of the presence of spinorial fields in the 4-dimensional 
membrane, we consider the Dirac massive equation. Since we are interested in massive local- 
ized states, we transform the equation of motion for fermions in a Shrodinger-like equation, a 
well-known procedure also used in our study of gravity localization in Sec. III. For fermions 
in a domain wall there is a similar analysis 37]] . Double walls were studied in 38]. 

With the transformations a^iz) 



a Ln (z)e 2Ap and a Rn (z) = a Rn (z)e 2Ap , Eqs. (E 
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FIG. 4: Plots of V p R {z) (right chiral mode potential) with (a) / = 0.1 (left) and (b) / = 2 (right) 
fixed. Curves correspond to p = 1 (tinner line), 3, 5 and 7 (thicker line). 



and (|25|) result in 

[-d 2 z + V p L }a Ln (z) = m 2 a Ln (z), (36) 
[-d 2 + V p R ]a Rn (z) = m 2 a Rn (z) , (37) 
where the Schrodinger potentials are 

V p L = -fd z cj> p e A » - f<j> p e A »d z A p + / 2 0> 2Ap (38) 
V p R = +fd z <P P e A - + f4> p e A "d z A p + f 2 <p 2 p e 2A ". (39) 

Due to the change on variables from y to z we have no explicit form for the potentials V p 
and V p . From the numerically known potential we can use the Numerov numeric method 
fis| to solve the Schrodinger equations. 

Fig. shows the potentials V p for p = 1, 3, 5 and fixed / = 0.1 and f — 2. Note from 
the figure that the potential for p — 1 is qualitatively different from p = 3, 5. For p = 1 and 
/ = 0.1 the potential has a maximum at the center of the brane, and the structure of the 
potential shows that there is no resonance. For p = 3,5, there appears a structure of two 
peaks separated by a local minimum at z — 0. The case where / = 2 is more interesting, 
where there appears a hole for p = 1. For p = 3, 5, 7 the height of the peaks is roughly 
the same, but the peaks are more distant for larger values of p, showing that the internal 
structure of the brane favor the appearance of resonances. 

Fig. (jSJ) shows the potential V p for several values of p. In this case there is only one 
minima for p — 1, where the potential is negative. This shows that a brane with this 
potential has a bound state, the zero-mode already studied. For larger values of p, this 
central minimum is separated into two minima located far from the brane for larger values 
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- 0, 1 o J 



FIG. 5: Plots of V p L (z) (left chiral mode potential) with (a) / = 0.1 (left) and (b) / = 2 (right) 
fixed. Curves correspond to p = 1 (tinner line), 3, 5 and 7 (thicker line). 

of p. We note that these minima are less pronounced for larger values of p, disfavoring the 
localization of the zero-modes. This agrees with Fig. [3^, where we see that the solution 
for the zero- mode ot-Ln{y) for p — 1 is characterized by a higher peak in comparison to 
what found for larger values of p. Fig. (jSj) also shows that the potential is characterized 
by two positive maxima. This shows that for energies between zero and this maxima one 
can investigate the presence of resonances. Note also that for p > 3 and for small values 
of /, say / = 0.1, with the increasing of p, the height of the two peaks from the potential 
decrease and are more apart one from the other (see Fig. ([5^)). This is a similar effect 
was also noted in the Sec. Ill of this paper, when studying the Shrodinger potential in 
the transverse-traceless sector of metric fluctuations. In this case we would expect that 
higher values of p are less effective to trap the resonant massive modes. However, such small 
values of / do not satisfy the inequality (|29|) . and zero- modes are non-normalizable. On the 
contrary, for larger values of /, for instance / = 2, the height of the two peaks are almost 
constant for p > 3 and are also more apart one from the other (see Fig. (JSb))- Such effect 



where such characteristics are 



where also found in studying gravity fluctuations in ref. 
compared to the phase transition of complete wetting in condensed matter systems. In our 
case this lead us to expect that higher values of p are more effective to trap the resonant 
massive modes, increasing their lifetimes. 

We can analyze more closely the influence of the 5-dimensional Yukawa coupling. As an 
specific example we chosen to fix p = 3 and varied the parameter / for right-chiral fermions. 
Fig. [6]shows that for / = 0.1, where inequality (|29|) is violated, the potential is qualitatively 
different for the righer values of / analyzed, assuming negative values far from the brane. 
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FIG. 6: Plots of V p R (z) with p = 3, for / = 0.1 (upper left), 0.5 (upper right) , 1.1 (lower left), 2 
(lower right). 

For the larger values of / presented in the figure, the position of the two maxima is roughly 
the same, independent of /. The height of the maxima increase faster with /, showing a 
tendency for more resonances to appear as the coupling / increases. 

Note that for all potentials in Figs, (jlj) and (j5J), the Shrodinger potential asymptotes to 
zero, and there is no gap in the spectra. The potentials for both chiralities show significant 
changes when the parameter p changes from 1 to larger values. The appearance of internal 
structure is the important point and can reveal other physical aspects related to fermion 
localization. In order to investigate the confining of the massive chiral modes, we must solve 
Eqs. f )36|) and fl37|) . It is easy to see that for massive modes that equations can be rewritten 
as 

QQ + a Rn = [d z + f<p p e A f} \-d z + f(p p e Ap }a Rn = m 2 a Rn , (40) 
Q + Qa Ln = [-d z + f4> p e Ap }[d z + f4> p e Ap ]a Ln = m 2 a Ln , 

corresponding to a supersymmetric quantum mechanics scenario. Written the equations of 
motion in this form, tachyonic modes are clearly forbidden. 
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FIG. 8: Plots of am for p = 1 (left), and (b) p = 9 (right, doted line), p = 11 (right, solid line). 
We fix / = 1 and m = 0.4. 



We numerically solve Eqs. (I5rjj) - (15T|) and plot the solutions for a Rn and ollti for some 
values of p and / in Figs. jQ and (JED, respectively. A method for solving such equations 
can be found in Refs. 18|, |28|. In each figure we show the solutions for p = 1 on the left 
representing kink models already considered in the literature 37|, |38( , whereas on the right 
we shown the solutions on the deformed walls. In these figures we considered higher values 
of p where the difference between the background treated here and the p — 1 is clearer. 
Note that the solutions show a transition region for both chiralities near the center of the 
brane. Far from the brane the solutions have the character of plane waves, signifying that 
the fermion escapes from the brane. 
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VI. FERMIONIC RESONANCES AND DIRAC FERMIONS 



The changing of variables that produced the Shrodinger equations presented in Eq. (1401) . 
lead us to adopt a quantum mechanical interpretation for ain and a^ n . One importance 
of studying resonances is connected to their information of the coupling of massive modes 
and the brane, illustrating how the mechanism of fermion trapping is being processed. In 
our case, we can interpret \Na±(z)\ 2 as the probability for finding the massive mode in the 
position z, with iV a normalization constant. In this way, calculating P(m) = |iVo7±(0)| 2 as 
a function of the mass m, we are able to detect resonant modes as large peaks in the plot 
P(m) versus m. 

First of all we investigate right chirality, where there is no zero-mode. Fig. Hh shows that 
for the case p — 1 and small / there is no local minima for the potential, and resonances are 
absent. The presence of such minima for p > 3 as depicted in Fig. Hb shows that resonances 
possibly exist for a wider range of /. In particular for / = 2, Fig. Hb shows that the 
maxima of the potentials for p = 3,5, 7... are nearly the same, as well as the local minima. 
However, the region in the potential near the local minima is broader for larger values of 
p. This signals that for a large value of /, larger values of p are more effective in trapping 
the fermionic massive modes in comparison to smaller values of p. This must be compared 
to the richer structure of the energy density of the branes for larger values of p, as noted in 
[hJ. We found resonance peaks for even parity wavefunctions and analyzed the influence 
of the parameter p and coupling constant /. 

Some results for right chiral fermions are shown in Fig. [9] for p = 3, 5 and 7 and three 
values of /. Fig. [10] shows in detail some thin resonant peaks corresponding to Fig. [9l 
with the corresponding values of m 2 identified in Table I. For p = 3 and / = 0.5 we see 
that there is no resonance, as the width at half maximum Am of the peak is bigger than 
the mass m corresponding to the peak; for / = 1.1 and / = 2 there appears one and three 
resonances, respectively. This agrees with our expectation that larger values of / favor the 
presence of resonances. Note that in general, for fixed p, larger values of / lead to larger 
number of peaks. Also, for fixed /, larger values of p corresponds to the thicker peaks with 
almost the same mass. The first peak is the thinnest, with corresponding larger lifetime. 
The increasing of the parameter p tends to increase the mass of the first resonance (this is 
more evident for larger values of /), whereas the increasing of / turns richer the spectrum, 
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FIG. 9: Plots of |P(0)| 2 versus m 2 with p = 3 (first line), p = 5 (second line) and p 
line). Coupling parameters are / = 0.5 (left figures), / = 1.1 (middle figures), and / 
figures), for even parity wavefunctions of fermions with right chirality. 

with more resonances. 



Right 


/ = 0.5 


/ = !•! 


/ = 2 


p = 3 


absent 


0.26061 


2.76451953; 5.2062; 6.94 


p = 5 


0.042 


0.0873678916; 0.4764 


4.6115419;5.8111; 6.746 


p=7 


0.0230 


0.26301844; 0.5134 


4.99069057;5.782825;6.476;6.99106 



TABLE I: First resonance peaks, even parity modes, right chirality. The table shows the corre- 
sponding values ofm 2 . 

Next we analyzed left chiral fermions. Fig. [11] shows the resonance peaks for p = 3, 5, 7 
with same values of / used previously, but now for left chirality. Fig. [12] shows in detail 
some thin resonant peaks corresponding to Fig. [TT], with the corresponding values of m 2 
identified in Table II. We found the same effects related to the influence of p and /. In 
particular, note that the mass corresponding to the broader resonance peak did not change 
significantly with p. One remarkable point observed with the increase with p is the tendency 
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= 7 (third 
= 2 (right 



FIG. 10: Detailed plots of some right chiral resonance peaks for |P(0)| 2 versus m 2 , corresponding 
to / = 1.1 (left) and / = 2 (right) from Fig. M 

for the resonances corresponding to tiny peaks to accumulate near the broader peak. 





/ = 0.5 


/ = 1.1 


/ = 2 


p = 3 


absent 


0.646 


4.63126; 6.593 


p = 5 


absent 


0.298164; 0.640 


5.52762; 6.552 


p = 7 


0.0432 


0.623; 0.40240 


5.591930; 6.318; 6.872 



TABLE II: First resonance peaks, even parity modes, left chirality. The table shows the correspond- 
ing values of m 2 . 

Now we resume our findings for the resonance spectrum for both chiralities. First of all 
note that from Eqs. (j4"U|) the spectrum of the right and left chiralities are related. Indeed, 
the spectrum starts with zero mode with right chirality and even parity. The two resonances 
with same m 2 that succeed in the spectrum are the first odd parity left chiral mode and the 
first even parity right chiral mode. Next we also have two resonances with same m 2 : the 
first even parity left chiral mode and the first odd parity right chiral mode. We successively 
have in the spectrum even and odd parity wavefunctions for left and right chiral modes with 
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FIG. 11: Plots of |P(0)| 2 versus m 2 with p = 3 (first line), p = 5 (second line) and p = 7 (third 
line). Coupling parameters are / = 0.5 (left figures), / = 1.1 (middle figures), and / = 2 (right 
figures), for even parity wavefunctions of fermions with left chirality. 



same values of m 2 . One can easily check the formation of Dirac fermions after studying 



the odd parity wave 
known procedure 



'unctions. This needs the changing the normalization procedure in a 
24J ] (for an analysis with models with two scalar fields, see Ref. 28]). We 
checked the correspondence of the spectrum for some values for resonances with odd parity 
wavefunctions, but we will not pursue this subject in detail, since the numerical procedure 
was tested with confidence and the formation of the Dirac fermions are guaranteed by the 
super symmetric quantum mechanics structure. Also the study of even parity wavefunctions 
for massive fermions is interesting since we can compare our findings with the results for 
massive gravitons in the first part of this paper. 

The spectra for p = 3, 5, 7 and / = 2 are depicted in Fig. [131 There one can see 
how the increasing in p increases the number of resonances. Note also that the resonances 
appear coupled in pairs for left and right chiralities (here we consider only even parity 
wavefunctions). We also see that an increasing of p increases the mass of the first resonance, 
generally being a very thin line in the spectrum. The fact that the lines for left and right 
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FIG. 12: Detailed plots of some left chiral resonance peaks for |P(0)| 2 versus m 2 , corresponding 
to / = 1.1 (left) and / = 2 (right) from Fig. ED 



3 
2 



FIG. 13: Spectrum for p = 3 (left), p = 5 (middle) and p = 7 (right) with / = 2 fixed, constructed 
with even parity modes. Black lines for left chiral massive modes. Grey lines are for right chiral 
ones. 

chiralities alternate each other in the spectrum is very important for capture the thinner 
lines in the numerical procedure. 
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VII. CONCLUSIONS 



In this work we have studied localization of fermions and gravity in a model of deformed 
branes. The parameter p of the model controls important features as the brane thickness 
and energy distribution along the extra dimension. We found that the increasing of p 
cause a splitting in the Ricci scalar, evidencing the appearance of an internal structure, 
as can be found also from the energy density analysis (see 10|). We investigate metric 
perturbations as decoupled from the scalar ones in the transverse-traceless gauge. From the 
asymptotic behavior of the Shrodinger-like equation it can be proved that gravity is localized 
for all parameters p, and tachyonic gravity modes are absent. Resonances in the gravity 
sector are observed for p = 1 and 3 as peaks in the \H^(0)\ 2 distribution. However, larger 
values of p lead to produce broader peaks that cannot technically characterize themselves as 
resonances. This means that thinner branes are more effective for trapping gravity. Next we 
studied the presence of fermionic zero-modes after introducing a Yukawa coupling between 
the 5-dimensional spinor and the scalar field, depending on a parameter /. After a chiral 
decomposition we obtained a Shrodinger-like equation also for fermionic scalar functions that 
characterize localization of left or right chiral modes and investigated its normalizability. We 
found that zero-mode exist only for left chirality and that tachyonic fermionic modes are 
absent. The investigation of fermionic resonances is a subject possible here only numerically. 
However, from the qualitative character of the Shrodinger potentials we were able to observe 
that for large values of / the behavior of localization of fermions is favored for larger values of 
p. This is contrary to the observed for gravity, where the localization of KK gravity modes are 
favored for small values of p. The analysis of small values of / could lead to similar behavior 
for gravitons and fermions. However, the parameter / cannot be reduced below a threshold 
where the presence of fermionic zero-modes are forbidden. For / above this threshold the 
effective action is finite. The presence of fermionic resonances was investigated for even 
parity wavefunctions both for left and right chiralities. We analyzed how the presence of 
the Yukawa coupling influenced our findings, with larger values of / and p favouring the 
increasing of resonances. Some resonances are extremely thin and more difficult to be found. 
The quantum mechanical supersymmetric character of the Shrodinger-like potentials for left 
and right chiral fermions guarantee that Dirac fermions are realized in the model and that 
the spectral lines for left and right fermions with same parity must alternate in the diagram. 



21 



This is very helpful in the numerical process for finding resonances, serving as a guide for 
reducing the numerical step when a resonance line is not captured. The lightest modes in 
the spectrum correspond to the thinnest lines, showing that those modes couple strongly 
with the brane in comparison to the KK modes with higher masses. This is expected as 
lightest modes have lower energy to escape from the brane Shrodinger potential. As the 
inverse of the peak width to half maximum is proportional to the lifetime of the resonance 



39(, thinner peaks may correspond to particles in nature with sufficiently large lifetimes to 



be important in phenomenology. 
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